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Abstract 


Radiative scattering by clouds leads to errors in the re- 
trieval of column densities and concentration profiles 
of atmospheric trace gas species from satellites. The 
SCanning Imaging Absorption SpectroMeter for Atmo- 
spheric CHartographY (SCIAMACHY) instrument on 
the Envisat satellite, principally designed to retrieve trace 
gases in the atmosphere, is also capable of detecting 
clouds. Cloud fraction is retrieved using the Optical 
Cloud Recognition Algorithm (OCRA). In this paper, this 
Level 2 operational cloud product is compared to the re- 
sults from other cloud algorithms such as FRESCO (Fast 
Retrieval Scheme for Clouds from the Oxygen A-band). 
It is shown that the retrieved cloud fractions correlate well 
between OCRA and FRESCO since the 28"" April 2004. 
FRESCO provides a consistent set of cloud products by 
retrieving simultaneously cloud fraction and cloud top 
pressure. Moreover, FRESCO retrievals are corroborated 
with cloud top pressure from MODIS. Other algorithms 
such as HICRU (cloud fraction) and SACURA (cloud 
top pressure) are also available to validate independently 
these cloud products. 
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1. INTRODUCTION 


SCIAMACHY (SCanning Imaging Absorption Spec- 
troMeter for Atmospheric CHartographY) is a space- 
borne spectrometer that flies on Envisat since March 
2002. The instrument is a national contribution by 
Germany-DLR, the Netherlands-NIVR and Belgium- 
BUSOC to the European Space Agency satellite En- 
visat. SCIAMACHY measures the solar radiation re- 
flected from the atmosphere in the wavelength range be- 
tween 240 and 2380 nm. This is recorded at relatively 
high resolution (0.2 to 1.5 nm) over the range 240 to 1750 
nm, and in selected regions between 2000 and 2400 nm. 
SCIAMACHY is extended as compared to its precursor 
GOME (Global Ozone Monitoring Experiment; [1]) with 
a wavelength range extending in the near IR region and 
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able to operate with new viewing geometries namely limb 
and sun, moon occultations. The primary scientific ob- 
jective of SCIAMACHY is the global measurement of 
various trace gases on a global scale [2]. 

The trace gases that SCIAMACHY detects, like O3, 
NO», BrO, SO, CO and CHy, occur not only in the 
stratosphere but also in the troposphere, where clouds 
reside. To detect these trace gases accurately, the 
presence and properties of clouds have to be known, 
in order to correct for the effect of clouds. The three 
cloud parameters that most strongly influence trace gas 
measurements are: cloud fraction, cloud albedo or cloud 
optical thickness, and cloud top pressure. Microphysical 
cloud parameters, like particle size and shape, are of 
minor importance [3]. 


Clouds affect the path of photons through the atmosphere 
and therefore change the interpretation of the depth of 
an absorption band. Taking ozone as an example, the 
effect of clouds on ozone retrieval can be regarded as 
due to two main effects [3]: (1) albedo effect: clouds 
act as a reflecting boundary below the ozone layer and 
enhance the depth of the ozone absorption bands, and (2) 
ghost column effect: clouds shield tropospheric ozone 
from observation. A third, smaller effect exists, which 
is the enhancement of the photon path inside clouds, 
causing an enhancement of the absorption line depth 
inside clouds [4]. 


In order to correct for the albedo effect of clouds, it is 
necessary to retrieve the cloud fraction and cloud albedo 
(or cloud optical thickness) or, equivalently, the product 
of cloud fraction and cloud albedo, which we call the 
effective cloud fraction. In order to correct for the ghost 
column effect of clouds, it is necessary to retrieve both 
the effective cloud fraction and the cloud top pressure. 


The above cloud effects not only occur for O3 but also 
for other gases. Especially for NO2, much of which is 
residing in the troposphere in polluted circumstances, the 
ghost column effect is very important [5]. Therefore, a 
cloud detection algorithm suitable for trace gas retrievals 
for SCIAMACHY should be able to detect at least 
effective cloud fraction and cloud pressure. If the cloud 
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pressure is available, then the tropospheric ozone column 
can be determined by comparing cloudy and clear pixels, 
using the convective-cloud-differential technique [6]. 


This paper is devoted to the SCIAMACHY cloud 
products’ validation. The operational Level 2 cloud 
fraction from OCRA is investigated and compared with 
the retrievals obtained from FRESCO. Moreover, the 
scientific cloud top pressure product from FRESCO 
is compared to results from the O2-O2 method and 
from another satellite instrument, namely MODIS on 
EOS/Terra. 


2. Cloud retrieval techniques for SCIAMACHY 


For reasons of timeliness and co-location, it is most prac- 
tical to retrieve cloud parameters from the SCIAMACHY 
data itself. To this aim various SCLIAMACHY channels 
and atmospheric absorption bands can be chosen. In the 
case of GOME, the Polarisation Measurement Devices 
(PMDs) had been used by several groups to detect cloud 
fraction. The advantage of the PMDs is that they have 
a better spatial resolution than the spectral channels. 
Reference [7] have developed various cloud algorithms 
for GOME, among which the Optical Cloud Recognition 
Algorithm (OCRA) technique. This technique uses the 
color of PMD images to select cloud free scenes for a 
global database of cloud free reflectances. The difference 
between the current PMD reflectances and this cloud 
free database is used to retrieve cloud fractions. To 
detect cloud pressure, however, the PMDs cannot be 
used because they have insufficient spectral resolution 
to detect absorption lines of well-mixed absorbers 
like oxygen. In the SCIAMACHY range there are 
several oxygen or oxygen-oxygen absorption bands. 
The strongest band is the O2 A-band at 760 nm. Fig- 
ure 1 shows the absorption lines of Oxygen in this region. 
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Fig. 1. The Oxygen molecular absorption optical thickness in the 
Oxygen A band of the Mid-Latitude Summer atmosphere. 


There is also the weaker O2 B-band around 680 nm, 
and several weak O2-O2 lines. The Initial Cloud Fitting 
Algorithm (ICFA), which is used in the operational 


data processing of GOME, uses the O2 A-band [8], [9]. 
ICFA only produces the cloud fraction; cloud pressure is 
taken from a climatology. Validation of ICFA has been 
performed by [9], and it was shown that ICFA produces 
inaccurate cloud fractions. The GOME-CAT (Cloud 
retrieval algorithm) method [10] uses a combination of 
PMD data and O, A-band data to retrieve cloud fraction, 
cloud optical thickness and cloud pressure. In the results 
of this algorithm, however, the cloud optical thickness 
appears to be unreliable. 


Recently, an intercomparison of different cloud retrieval 
algorithms that exist for GOME was published by [11]. 
They compared ICFA, OCRA and FRESCO cloud 
fractions with groundbased synoptical observations. 
It was shown that the effective cloud fraction as used 
by FRESCO underestimates the real cloud fraction for 
optically thin clouds, as expected, but correlates well 
with synoptical cloud observations. Other recent cloud 
algorithm comparisons for GOME and SCIAMACHY 
have been performed. Reference [12] developed an 
advanced PMD technique for cloud fraction retrieval, 
called HICRU (Heidelberg Iterative Cloud Retrieval 
Utilities). SACURA (SemiAnalytical CloUd Retrieval 
Algorithm) from [13] determines the cloud top pressure 
and geometrical thickness using measurements of the 
cloud reflection function. 


Apart from the O2 A-band, also the 477-nm absorption 
band of the collision-complex O2-O2 has been used 
for retrieval of cloud fraction and cloud pressure from 
GOME [14]. The O2-O2 method will be used as opera- 
tional cloud detection method for OMI, to be launched 
on board of EOS-Aura in July 2004. The O2-O2 method 
is more sensitive to low clouds than to high clouds, due 
to the pressure-squared dependence. 


Another atmospheric scattering process that can be used 
for cloud pressure retrieval is Raman scattering. Raman 
scattering is the inelastic part of Rayleigh scattering by 
nitrogen and oxygen, which fills-in solar Fraunhofer 
lines in scattered light from the Earth. The amount of 
filling-in is more or less proportional to the pressure 
of the cloud. Recently, this Raman method has been 
successfully applied to GOME data [15]. 


In cloud retrievals from satellite, it is very important to 
have a good estimate of the surface albedo. The reason is 
that cloud detection is usually performed by comparing 
the measured reflectivity with the expected reflectivity 
from a cloud-free scene. Recently, a new spectral surface 
albedo database based on 5.5 years of GOME data has 
been published [16]. This database is freely available 
and is used in FRESCO in this study. 


2.1 FRESCO method 


The FRESCO (Fast Retrieval Scheme for Cloud Ob- 
servables) method has been developed by [17] as a 
simple, fast, and robust algorithm to provide cloud 


information for cloud correction of ozone. FRESCO 
uses the reflectivity in three 1-nm wide windows of the 
O2 A-band: 758-759 nm, 760-761 nm, and 765-766 
nm. The measured reflectivity is compared to a modeled 
reflectivity, as computed for a simple cloud model. In 
this model the cloud is assumed to be a Lambertian 
reflector with albedo 0.8 below a clear atmosphere, in 
which only O2 absorption is taken into account. To 
simulate the spectrum of a partly cloudy pixel, a simple 
atmospheric transmission model is used, in which the 
atmosphere above the ground surface (for the clear 
part of the pixel) or cloud (for the cloudy part of the 
pixel) is treated as a purely absorbing, non-scattering, 
medium. The retrieved parameters are the effective 
cloud fraction (between O and 1) and the cloud pres- 
sure. Absorption by oxygen inside the cloud is neglected. 


FRESCO is used in the fast-delivery processing of 
GOME ozone data [18] and in the SCIAMACHY ozone 
processor TOSOMI, which provides ozone data from 
SCIAMACHY within the ESA TEMIS project (see 
http://www.temis.nl). 


FRESCO has been validated regionally and globally 
[17], [13], [11], and is proven to be reliable. However, 
FRESCO overestimates the cloud fraction retrievals 
over the Sahara region. This is under investigation but a 
possible explanation could be the higher albedo over this 
desert, from sand and aerosols, than in the low resolution 
surface albedo database used by FRESCO. Secondly, 
FRESCO overestimates cloud top pressure by on average 
about 50 hPa, as compared to infrared techniques. A 
possible reason is the fact that photon paths inside the 
cloud enhance the O2 A band depth. We note that 
for cloud correction of ozone and other trace gases, a 
qualitatively similar path enhancement is expected. 


2.2 OCRA method 


OCRA was developed by [19]. It is a PMD algorithm 
using the ”red”, ” green” and blue” PMDs for cloud frac- 
tion determination. The cloud fraction is calculated from 
a linear interpolation of the current PMD reflectance 
between a minimum (cloud free) value and a maximum 
(totally cloudy) value, the latter being fixed by scaling 
factors. An average of the results of the three PMDs is 
used. Reflectances of cloud free scenes are stored on a 
global grid. They are determined once from a large set 
of historic data using color distances between measured 
PMD color and the white point in the chromaticity 
diagram. Reflectances with the largest distance to the 
white point are selected to represent the cloud free case 
at a certain grid point. Currently, a database calculated 
from GOME PMD measurements is applied. For the 
totally cloudy case, one scaling factor is used for each 
PMD, representing the maximum reflectance. 


3. VALIDATION RESULTS 


In this section we directly intercompare the retrieved 
cloud paramaters from the different cloud algorithms. 
More specifically, results obtained by applying the 
FRESCO algorithm to SCIAMACHY data are compared 
with the equivalent Level 2 operational cloud fraction 
from OCRA. Moreover, the retrieved cloud top pressure 
results are assessed against another satellite instrument 
based values (MODIS). 


3.1 Cloud fraction 


Firstly, Figure 2 illustrates a comparison between OCRA 
and FRESCO cloud fraction retrievals for the validation 
reference dataset. This dataset corresponds to the states 
shown in Figure 3 covering the period from August to 
October 2002. A good agreement is obtained between 
OCRA and FRESCO with a standard deviation of 0.1. 
The discrepancies at the cloud fraction value of 1.0 are 
due to FRESCO algorithm assumptions at high latitudes. 
Indeed, in these regions, the algorithm switches to 
snow/ice mode where an effective cloud fraction of 1.0 
is assumed. 
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Fig. 2. Comparison of cloud fraction derived from SCIAMACHY 
Data between OCRA (Scia) and FRESCO (Tosomi) retrieval cloud 
algorithms for the validation reference dataset. 


Secondly, a more specific comparison between OCRA 
and FRESCO cloud fraction retrievals is shown in 
Figure 4. This concerns the SCIAMACHY data from the 
reference orbit 2510 (August 23, 2002). It is clear that 
both cloud fractions correlate well. The linear fit gives 
a determination coefficient (R?) of 0.95 with a slope of 
0.94. 


However, plotting the cloud fraction difference between 
OCRA and FRESCO as a function of the viewing zenith 
angle (Figure 5) for the same orbit produces a rather 
surprising parabolic response. It was found that this 
dependency in scanning angle is due to an inadequate 
definition of reflectances in OCRA, which takes the 
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Fig. 3. SCIAMACHY measurements locations corresponding to the 
validation reference dataset used in the comparison between OCRA and 
FRESCO cloud fractions of Figure 2. The measurements date from the 
period from August to October 2002. 
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Fig. 4. Comparison of cloud fractions derived from the OCRA and 
FRESCO algorithms using SCIAMACHY data of orbit 2510 (August 
23, 2002). The dotted line is the one-to-one agreement while the solid 
blue line is the best fit produced by a regression analysis. 


cosine of this angle into account. 


Furthermore, in Figure 6 the retrieved cloud fractions 
from both FRESCO and the O2-O2 method are compared 
with OCRA values in the case of the orbit 2510, corrob- 
orating the previous results. 


The good correlation in the retrieved cloud fraction 
values between OCRA and FRESCO with the veri- 
fication orbit 2510, but also 2509, is also found for 
SCIAMACHY data on Near-Real-Time orbits since 
April 28, 2004 (Figures 7 and 8). Indeed, beforehand, 
an error with a Near-Real-Time Level 1-2 processor 
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Fig. 5. Cloud fractions difference between OCRA and FRESCO as 
a function of the scanning angle. This is derived from orbit 2510 of 
SCIAMACHY data (August 23, 2002). 
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Fig. 6. Comparison of cloud fraction derived from FRESCO (blue) 
and the O2-O2 method (red) against OCRA using SCIAMACHY data 
of orbit 2510 (August 23, 2002). 


initialisation file prevented OCRA to retrieve accurate 
values of cloud fraction. 


Another promising cloud algorithm (HICRU) is currently 
applied on SCIAMACHY data to retrieve cloud fraction. 
First results (Figure 9) showed that the retrieved cloud 
fractions agree well with the values obtained from 
GOME data and, also, with images from Meteosat [20]. 


3.2 Cloud top pressure 


The comparison between SCIAMACHY operational 
Level 2 and FRESCO cloud top pressures from the 
reference orbit 2510 (August 23, 2002) underlines a 
strong disagreement. Indeed, while FRESCO provides 
a consistent set of cloud parameters (simultaneous 
retrieval), the operational product is derived from the 
ISCCP monthly mean cloud top height database, which 
is a climatology based on measurements of clouds 
from satellites (MeteoSat, GMS, GOES and NOAA). 
However, in a future version of the operational processor, 
an oxygen band algorithm for cloud top pressure retrieval 
will be included. 
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Fig. 7. Comparison of cloud fraction derived from the OCRA and 
FRESCO algorithms using SCIAMACHY data of orbit 11326 (April 
30, 2004). The dotted line is the one-to-one agreement while the solid 
blue line is the best fit produced by a regression analysis. 
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Fig. 8. Comparison of cloud fraction derived from the OCRA and 
FRESCO algorithms using SCIAMACHY data of orbit 11568 (May 17, 
2004). The dotted line is the one-to-one agreement while the solid blue 
line is the best fit produced by a regression analysis. 


This is corroborated in the comparison of the SCIA- 
MACHY operational product (Figure 10) and FRESCO 


Fig. 9. Comparison of cloud fraction derived from the HICRU algo- 
rithm using SCIAMACHY (left) and GOME (right) data from January 
24, 2003. 


(Figure 11) retrieved results with MODIS [21] co-located 
values for the orbit 2510. MODIS Level 2 granules over- 
lapping the SCIAMACHY 2510 orbit are processed. 
While the operational product exhibits important differ- 
ences with MODIS values, FRESCO shows a reasonable 
agreement with differences about 100 hPa. 
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Fig. 10. Comparison of the SCIAMACHY operational cloud top 
pressure with MODIS co-located values for the orbit 2510 (August 23, 
2002). The dotted line is the one-to-one agreement while the solid line 
is the best fit produced by a regression analysis. 


Furthermore, in Figure 12, the retrieved cloud top 
pressures from FRESCO (blue), the O2-O2 method 
(red) and the operational product (black) are com- 
pared with MODIS co-located values in the case of 
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Fig. 11. Comparison of cloud top pressure derived from FRESCO 
using SCIAMACHY data with MODIS co-located values for the orbit 
2510 (August 23, 2002). The dotted line is the one-to-one agreement 
while the solid line is the best fit produced by a regression analysis. 


the orbit 2510. For FRESCO and the O2-O2 method, 
only the retrieved values concerning optically thick 
clouds are considered (cloud optical thickness of 5 
or more). This improves the correlation of FRESCO 
with MODIS obtained in Figure 11 and confirms the 
slight overestimation in FRESCO retrieved values. To 
complete this validation, there are also other algo- 
rithms which should be employed such as SACURA 
(see http://www.ww.iup.physik.uni-bremen.de/scia-arc) 
providing cloud top pressure retrieval (dataset not yet 
available on-line) using SCIAMACHY data [13]. 


4. CONCLUSION 


To validate SCIAMACHY cloud products, an intercom- 
parison has been undertaken between the results of the 
operational algorithm OCRA and other algorithms such 
as FRESCO or the O2-O2 method. Moreover, the val- 
idation has been extended with MODIS co-located data 
in the case of the cloud top pressure. This shows that 
OCRA operational Level 2 cloud fraction product cor- 
relates well with FRESCO and the O2-O2 method us- 
ing SCIAMACHY data. The accuracy of this product 
is about 0.1, as compared to FRESCO, and could even 
be improved by correcting from the scanning angle de- 
pendency. However, for the current implementation of 
the OCRA method in the operational SCIAMACHY pro- 
cessor, this is the first, and encouraging validation study. 
Further studies will be needed to define OCRA opti- 
mised parameter settings. In terms of cloud top pressure, 
FRESCO matches well the MODIS values with an ac- 
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Fig. 12. Comparison of cloud top pressure derived from FRESCO 
(blue), the O2-O2 method (red) and SCIAMACHY operational product 
(black) against the MODIS co-located values for the orbit 2510 (August 
23, 2002). 


curacy about 100 hPa while SCIAMACHY operational 
product is derived from a climatological database. There- 
fore, FRESCO provides a consistent set of cloud param- 
eters by its simultaneous retrieval. The validation of 
these cloud products should be completed by using other 
available algorithms such as HICRU (cloud fraction) or 
SACURA (cloud top pressure) which evaluate indepen- 
dently cloud fraction and cloud top pressure. 
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